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Abstract 

We present a fluid dynamics video showing the results of a Q-bilUon 

atom molecular dynamics simulation of complex fluid flow in molten 
copper and aluminum. Starting with an atomically flat interface, a 
shear is imposed along the copper-aluminum interface and random 
atomic fluctuations seed the formation of vortices. These vortices grow 
due to the Kelvin-Helmholtz instability. The resulting vortical struc- 
tures are beautifully intricate, decorated with secondary instabilities 
and complex mixing phenomena. 

This work performed under the auspices of the U.S. Department of 
Energy by Lawrence Livermore National Laboratory under Contract 
DE-AC52-07NA27344. 

1 Introduction 

In nature, fluids are comprised of atoms and molecules. The complex struc- 
tures that form in fluid flow arc the result of the collective motions of those 
elementary constituents. In simulation, computers struggle to deal with the 
multitude of atoms needed to represent complex fluid flows, but as comput- 
ers become faster and more capable it is now possible to perform molecular 
dynamics simulations sufficiently large to reproduce fluid dynamics phenom- 
ena. 
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Our video, (large and small) shows the formation and evolution of vor- 
tices formed by shear flow in a 9-billion atom molecular dynamics simulation. 
The simulation domain is a quasi-2D geometry 12|Um x 6//m x 2nm with an 
equal number of aluminum and copper atoms at a temperature of 2000 K, 
well above the melting temperature. Energies and forces between atoms are 
calculated using the EAM potential[l] of Mason, Rudd, and Sutton[2]. The 
simulation was performed on the BG/L computer at LLNL using 212,000 
CPUs. The calculation required over 36 x 10^ CPU hours or roughly 4 CPU 
millennia. Further details of the computational aspects of this work can be 
found in ref [3J. 

The simulation begins with laminar flow in which the copper at the 
bottom of the image flows uniformly to the left and the aluminum at the 
top flows to the right with an atomically flat interface between the two. The 
relative velocity of the two layers is 2000 m/s. Color in the video indicates 
the fraction of copper atoms (red is 100% Cu; purple is 100% Al). 

The simulation runs for 10^ time steps or a total time of 2 ns during 
which vortices develop and grow at the interface due to the Kelvin-Helmholtz 
instability. As time progresses we see the vortices combine and grow as well 
as the formation of secondary instabilities. Work to compare the results of 
this simulation to traditional hydrodynamic methods is ongoing. 
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